Quantum simulator for the Ising model with electrons floating on a helium film 
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We propose a physical setup that can be used to simulate the quantum dynamics of the Ising 
model in a transverse field with present-day technology. Our scheme consists of electrons floating on 
superfluid helium which interact via Coulomb forces. In the limit of low temperatures, the system 
will stay near its ground state where its Hamiltonian is equivalent to the Ising model and thus shows 
phenomena such as quantum criticality. Furthermore, the proposed design could be generalized in 
order to study interacting field theories (e.g., A^ 4 ) and adiabatic quantum computers. 



PACS numbers: 03.67.Ac; 75.10.Hk. 

Introduction Richard Feynman's observation [1] that 
classical computers cannot effectively simulate quantum 
systems bred widespread interest in quantum computa- 
tion. He thought up the idea of a quantum processor 
which uses the effects of quantum theory instead of classi- 
cal physics. As an example, Feynman proposed a univer- 
sal quantum simulator consisting of a lattice of spins with 
nearest neighbor interactions that are freely specifiable 
and can efficiently reproduce the dynamics of any other 
many-particle quantum system with a finite-dimensional 
state space 1]. Although such universal quantum com- 
puters of sufficient size (e.g., number of QuBits, i.e., 
spins) are not available yet, it is possible to design a spe- 
cial quantum system in the laboratory which simulates 
the quantum dynamics of a particular model of interest. 
Such a designed quantum system can then be regarded 
as a special quantum computer (instead of a universal 
one, which is more challenging) which just performs the 
desired quantum simulation, see, e.g., [1, SHI- 

In the following, we present a design for a quantum 
simulator for the Ising spin chain in a transverse field 
and demonstrate that it could be feasible with present- 
day technology, i.e., electrons floating on a thin super- 
fluid Helium film. A similar idea based on trapped ions 
has been pursued in Q. Nevertheless, since different ex- 
perimental realizations possess distinct advantages and 
drawbacks, it is still worthwhile to study an alternative 
set-up. For example, the number of coherently controlled 
ions in a trap is rather limited at present, whereas our 
proposal can be scaled up to a large number of electrons 
more easily - which is important for exploring the con- 
tinuum limit and scaling properties etc. 

The model We want to simulate the quantum dynam- 
ics of the one-dimensional Ising chain consisting of n 
spins with nearest-neighbor interaction J plus a trans- 
verse field r along the x-direction (h = 1) 



H 



E 



(1) 



where a-j = (a* , ctJ , <rj) are the spin-1/2 Pauli matrices 
acting on the jth qubit. This model has been employed in 



the study of quantum phase transitions and percolation 
theory JH] , spin glasses [f| , as well as quantum anneal- 
ing [7|, LSI] etc. Although the Hamiltonian ([I]) is quite sim- 
ple and can be diagonalized analytically, the Ising model 
is considered a paradigmatic example [5( for second-order 
quantum phase transitions and is rich enough to dis- 
play most of the basic phenomena near quantum critical 
points. For T ^> J, the ground state is paramagnetic 
| — >—>■—>■ . . . ) with all spins polarized along the x axis. 
In the opposite limit r <C J, the nature of the ground 
state(s) changes qualitatively and there are two degen- 
erate ferromagnetic phases with all spins pointing either 
up or down along the z axis | tTT • ■ • ) or I Hi • ■ ■ )■ The 
two regimes are separated by a quantum phase transi- 
tion at the critical point r cr = J, where the excitation 
gap vanishes (in the thermodynamic limit n f oo) and 
the response time diverges. As a result, driving the sys- 
tem through its quantum critical point at a finite sweep 
rate entails interesting non-equilibrium phenomena such 
as the creation of topological defects, i.e., kinks [iij]. Fur- 
thermore, the transverse Ising model can also be used 
to study the order-disorder transitions at zero temper- 
ature driven by quantum fluctuations [B|, 0]. Finally, 
two-dimensional generalizations of the Ising model can 
be mapped onto certain adiabatic quantum algorithms 
(see, e.g., [l3|). However, due to the evanescent excita- 
tion energies, such a phase transition is rather vulnerable 
to decoherence, which must be taken into account [ill ]. 

The analogue In order to reproduce the quantum dy- 
namics of the 1+1 dimensional Ising model (p}, we pro- 
pose trapping a large number of electrons on a low- 
temperature helium film of thickness h (e.g., ft = 110 
nm) adsorbed on a silicon substrate [l2| ■ Due to the po- 
larizability e ss 1.06 of the Helium film, the electrons are 
bound to its surface (i.e., in z-direction) via their image 
charges and the large potential barrier (around 1 eV) for 
penetration into the helium film [13]. Since the binding 
energy of around 8 K is much larger than the tempera- 
ture T (below 1 K) and the width of the electron wave 
packet in z-direction (of order 8 nm) is much smaller 
than all other relevant length scales, the electron motion 
is approximately two-dimensional (x, ?/-plane). 
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FIG. 1: Sketch of the proposed analogue quantum simulator. 
Electrons (e~) are floating on a low-temperature helium film 
of height h adsorbed on a silicon substrate. A double-well 
potential for each single electron is created by a pair of golden 
spheres of radius a and distance d on the bottom of the helium 
film. The double wells at each site provide two lowest states 
of the electron and model the spin states | j) and | J.) at each 
site j. The tunneling rate between the two wells corresponds 
to the transverse field term Tcrf. The electrons are lined up 
at distances A and interact via Coulomb forces, which creates 
the term J<r|cr| +1 . 



In our scheme, each single electron on top of the he- 
lium film is trapped by a pair of golden spheres of radius 
a (e.g., a = 10 nm) and distance d (e.g., d — 60 nm) at- 
tached to the silicon substrate (i.e., on the bottom of the 
helium him, cf. Fig. [1]). Depending on its position x,y, 
the electron will also induce image charges in the two 
golden spheres (which act as a pair of quantum dots) 
and hence experience a double-well potential 



U w (x,y) 



ae 2 (x 2 +y 2 + a 2 + /3 2 ) /4 



TVS 



(x 2 +y 2 + a 2 + f3 2 ) - Aa 2 y 2 



(2) 



with a = d/2 + a and f3 2 = h 2 



Since this po- 



tential is quite deep and symmetric U(x,y) — U(x,—y), 
cf. Fig. [21 the ground state wave- function tps(x, y) is given 
by the symmetric superposition of the two Wannier states 
tpo(x, ±y) while the hrst excited state i/ja{x, y) is the anti- 
symmetric combination 



ips(x,y) 
iI>a(x,v) : 



. Mx,y) + Mx,-y) ^ I T) + II) 

i>o(x,y) - ip (x,-y) _^ | T) - 1 1) 
V2 V2 



(3) 



For a sufficiently high potential barrier between the two 
wells, the Wannier state tpo(x, y) is strongly concentrated 
in the left well and models the spin state | f) and vice 
versa. The tunneling between the two states is then de- 
scribed by the Pauli operator a x with a x \ f) = | J.) and 
° x \ i) = I t) such, that the tunneling rate, given by the 
difference of the eigenenergies Ea — Es of V'S and -0a, 
corresponds to the transverse field T in Eq. (JXJ) . In the 



FIG. 2: Sketch of the double-well potential U(y) with four 
turning points for the energy Eq. 



limit of strong localization (i.e., weak tunneling), the en- 
ergy splitting Ea — Es between the two levels can be 
estimated via the WKB approximation 



14] 



E\ — Es ~ — exp 

7T 



Wo 



dy \p(y)\ 



(10 



(4) 



Here u> is the oscillation frequency (within one well) and 
±yo are the two inner (classical) turning points, cf. Fig. [21 
The integrand is given by p(x, y) — yj2m e [E — U(x,y)], 
where we can set x = since the tunneling probability 
away from the x = 0-axis is strongly suppressed. Finally, 
the energy Eq determines the turning points and m e is 
the electron mass. For the parameters above, each valley 
can well be approximated by a harmonic oscillator 



U w (x,y w ±t/ min ) 



4tt£/3' 



.(X 2 + [y T y min ] 2 ) j ( 5 ) 



and thus we obtain Eq ~ \J ae 2 /2it£m e (3 4 rs uj. 

So far, we derived the term TctJ in Eq. {T]) via Eqs. ([2]), 
(|3|), (J4j), and ([5]). In order to simulate the remaining part, 
we propose to line up the pairs of quantum dots at equal 
distances A (e.g., A = 600 nm), where the parameters are 
supposed to obey the following hierarchy 



A> h > d> a . 



(6) 



In this limit, the interaction between the electrons will 
be dominated by the direct Coulomb repulsion between 
nearest neighbors U c (x,y) = X)J=i £^c" J+1 with n denot- 
ing the number of electrons floating on the helium film. 
For A 3> d, we may Taylor expand the Coulomb interac- 
tion into powers of y/X due to y ~ ±d/2. The zeroth- 
order term is constant and thus irrelevant while the first- 
order contributions vanish (up to boundary terms) after 
the sum over sites j. Thus, the leading term is bilinear 
in the electron positions 



U c (x,y) w — — 3^»i , (7) 

27T£ (A + a + 4a) ~{ 

and precisely corresponds to the Jcr^a Z j +1 term in Eq. ([T]) 
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with the effective coupling 

f e\d + 2af 

8^£ (A + d + 4a) 3 ' { > 

Experimental parameters For the example values 
given in the text, we obtain r ~ 0.1 K for the tunneling 
rate and the same value J w 0.1 K for the effective cou- 
pling, i.e., we are precisely in the quantum critical regime. 
However, deviations from this critical point should be 
easy to realize experimentally by varying the height h of 
the helium film, since the tunneling rate depends strongly 
(in fact, exponentially) on h, whereas the Coulomb force 
remains approximately constant. In order to see quan- 
tum critical behavior, i.e., to avoid thermal fluctuations, 
the temperature should ideally be well below this value 

0. 1 K (or at least not far above it). 

Furthermore, the Coulomb repulsion energy between 
two electrons (zeroth-order term) of about 11 K would 
tend to destabilize the electron chain. Fortunately, this 
effect is compensated by the binding energy between the 
electron and its image on the sphere, which is around 
13 K and thus stabilizes the electron chain. The proba- 
bility for the electron to penetrate the helium film by tun- 
neling to one of the golden spheres is extremely small (of 
order 10 -16 ) and can be neglected. Finally, the ground- 
state energy Eq w 1.4 K (within the harmonic oscilla- 
tor approximation) is reasonably well below the barrier 
height Uq w 3.1 K such that the WKB approximation 
should provide a reasonable estimate. (The tunneling 
probability of 0.08 is also small enough.) On the other 
hand, Eq s» 1.4 K is a measure of the distance between 
the two lowest-lying states in Eq. ([3]) and the remaining 
excited states in the double- well potential. As a result, 
these additional states do not play a role for tempera- 
tures well below one Kelvin and thus the Hamiltonian 
(fTj) provides the correct low-temperature description. 

Read-out scheme Having successfully simulated the 
Ising Hamiltonian fl} , one is lead to the question of how 
to actually measure its properties, e.g., how to detect 
signatures of quantum critical behavior. As one possi- 
bility, let us imagine enclosing the Ising chain symmetri- 
cally by two electrodes in the form of spheres of radius 
R = 100 /urn and a distance of 1 mm aligned along the 
chain axis. Applying a voltage of 1/xV, an approximately 
constant electric field of 4 x 10~ 4 V/m acts on all the 
electrons and induces the perturbation Hamiltonian 

n 

H peIt = ^ 7 o-?, (9) 

3=1 

corresponding to a longitudinal field (in addition to the 
transversal one r crj). For d = 60 nm, we get 7 s=s 0.1 fiK, 

1. e., a very weak perturbation 7 -C T. 

Deep in the paramagnetic phase T ^S> J, the response 
of the system to this weak perturbation 7 <C T is rather 
small (cr?) ~ 7/r. Approaching the phase transition, 



however, the static susceptibility x-y = lim 7 ^o(cJ)/7 
grows and finally diverges at the critical point. In the 
broken symmetry phase, the perturbation (J9J) lifts the 
degeneracy cr| — > — crj and hence the response is non- 
analytic, i.e., independent of the smallness of 7: e.g., 
for J 3> r, we have (crj) = sign(7) = ±1. This sig- 
nal (cr|) indicating the phase transition can be picked 
up by the two electrodes for which the Ising chain acts 
like a dielectric medium and induces a voltage shift of 
order nano-Volt per electron (for J>r), which should 
be measurable for a sufficiently large number of sites. In 
addition to the static case, one could also study the time- 
resolved response (crj (t)} to a varying voltage 7 (/;'), which 
is determined by the dynamical correlator (cr? (t')o'j (t)) 
in lowest-order response theory. 

Even in the absence of an externally imposed voltage, 
the chain induces spontaneous voltage fluctuations in the 
electrodes, which are strongest (of order nano-Volt per 
electron) deep in the ferromagnetic phase. The variance 
of these fluctuations yields the correlator sum X]ij( CT f a j) 
which is an order parameter for the phase transition and 
allows us to detect topological defects (i.e., kinks) which 
might have been produced during the sweep to the fer- 
romagnetic phase: In the presence of a kink, the ground- 
state signal J2ij(o~fo-j) — n 2 is drastically reduced (in 
average to Xaj(°f cr J) = "- 2 /3) depending on the kink 
position. If the kink is precisely in the middle of the 
Ising chain, we get a vanishing signal J2ij( a i a j) = 0' 
whereas a kink near the boundaries does not diminish 
the signal strongly. 

Disorder and decoherence In a real experimental set- 
up, the Hamiltonian will not be exactly equivalent to ([I]) 
due to imperfections such as electric stray fields, varia- 
tions in the film thickness h and further geometric param- 
eters a, d, and A etc. Therefore, the original expression 
^ will typically be altered to 

n 

a y. to °j + + tw> , (10) 

3=1 

where F, = F + £Tj and Jj = J + 5Jj. Assuming that the 
disorder parameters STj, SJj, and jj are much smaller 
than the excitation gap A = 2\ J — T\ of the undisturbed 
system (in the continuum limit), the impact of these im- 
perfections will be suppressed. Near the critical point 
J « r, however, this argument fails. Still, for a fi- 
nite number n of electrons, one retains a minimum gap 
(within the symmetric or anti-symmetric subspace, re- 
spectively) of order J/n. Exploiting this gap might be 
suitable for a reasonably small systems, but for n > 100 
electrons, the required accuracy on the sub-percent level 
is probably hard to achieve experimentally. E.g., decreas- 
ing the diameter of the golden spheres by ten percent 
with the other values remaining the same as before, the 
tunneling rate increases by fifty percent. 

For a sufficiently large number of electrons, the disor- 
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der induced by imperfections will become relevant near 
the critical point (in one spatial dimension) in view of 
the critical exponent v = 1 of the Ising model, see, e.g., 
0. (I.e., the renormalization flow is directed away from 
the homogeneous situation.) In this case, one would ex- 
pect effects such as local paramagnetic regions inside the 
global ferromagnetic phase and percolation transitions 
etc. Therefore, turning this drawback into an advan- 
tage, one might generate these imperfections on purpose 
in order to study the impact of disorder onto the phase 
transition. In contrast to the original Hamiltonian |T]), 
the above form (fTO]) is no longer analytically solvable and 
hence much less is known about its properties. Finally, in 
a real set-up, the system will also experience decoherence 



due to the inevitable coupling to the environment 



These effects could be incorporated by operator-valued 
variations STj, SJj, and jj associated to the degrees of 
freedom of the environment - where the same arguments 
apply as before. 

Summary We have proposed a design for the simula- 
tion of the quantum Ising model with a system of elec- 
trons floating on a liquid helium film adsorbed on a silicon 
substrate. Since the energy level splitting (tunneling rate 
r) depends exponentially on the thickness of the helium 
film h, we may tune the system through the quantum 
phase transition by changing h - which might even be 
feasible in a time-dependent manner, cf. Q. The created 
topological defects (kinks) could be detected via a strong 
reduction of the spontaneous voltage fluctuations in com- 
parison with the homogeneous ferromagnetic phase. 

Furthermore, a suitable generalization to two spatial 
dimensions might be relevant for adiabatic quantum al- 
gorithms, see, e.g., (loj . Note that the realization of a 
sequential quantum computer based on a set of electrons 
floating on a helium film has been proposed in [l5j ]. In 
contrast, our proposal is not suited for universal compu- 
tations, but (as one would expect) should be easier to 
realize experimentally. 

Exploring a different limit, where many eigenstates of 
the double-well potential contribute, the proposed set- 



up could simulate the lattice version of interacting field 
theories such as the A0 4 -model in 1+1 dimensions. 
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